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We investigate the spatial motion of the trapped atom with the electromagnetically induced trans- 
parency (EIT) configuration where the two Rabi transitions are coupled to two classical light fields 
respectively with the same detuning. When the internal degrees of freedom can be decoupled adi- 
abatically from the spatial motion of the center of mass via the Born-Oppenheimer approximation, 
it is demonstrated that the lights of certain profile can provide the atom with an effective field of 
magnetic monopole, which is the so-called induced gauge field relevant to the Berry's phase. Such 
an artificial magnetic monopole structure manifests itself in the characterizing energy spectrum. 



PACS number:03.65.Vf,42.50.Vk,14.80.Hv 

Modern concept of magnetic monopole in quantum me- 
chanics was postulated by Dirac in 1931 [1]. Since that 
time physicist have been making efforts to seek for the 
magnetic monopole in real space for more than seventy 
years. Though convincing evidence for its existence in 
real space has not yet been found, the theoretical concep- 
tion of magnetic monopole has initiated many important 
progresses in both physics and mathematics. In fact, in 
the extremely high energy scale that we have not reached 
at present, the grand unification theory [2] predicted the 
magnetic monopole as a consequence of the beautiful 
topology structure of Yang-Mills theory [3]. The dis- 
covery of Berry phase also resulted in a physical imple- 
mentation of magnetic monopole in the parameter space 
[4]. Precisely speaking, the degeneracy point in the pa- 
rameter space acts like a magnetic monopole caused by 
an effective gauge field. The Berry phase based magnetic 
monopole of this kind can be demonstrated in association 
with the anomalous Hall effect of ferromagnetic metals 
[5] where the slowly changing parameter is just the crys- 
tal momentum. Now we consider an artificial realization 
of magnetic monopole in real space. 

For the quantum adiabatic process induced Berry 
phase, it is well known that when the slowly varying pa- 
rameters are the dynamic variables of a subsystem inter- 
acting with another subsystem with fast varying vari- 
ables, the adiabatic separation of the two subsystems 
via the Born-Oppenheimer approximation [7] can provide 
the slow subspace with a scaler and a vector potentials 
[6] called the induced gauge potential. In the neutron 
spin precession experiment [8] the Aharonov-Borhm effect 
caused by this vector potential was pointed out as a man- 
ifestation of Berry phase. In this article, we will derive 
the monopole type induced gauge field for the spatial mo- 
tion of A— type atom interacting with control and probe 
laser beams which drive the transitions |e)-|l) and|e)-|2) 
respectively (see Fig. 1) [9,10]. We will also show that 
when the atom is cold enough, its internal degrees of 
freedom can decouple adiabatically from the spatial mo- 
tion of the center of mass. Correspondingly, the Born- 
Oppenheimer approximation provides the atomic center 



of mass with an effective magnetic monopole field as a 
special induced gauge field relevant to the Berry phase 
when the lights are artificially shaped in certain profiles. 
We predict that such an artificial magnetic monopole can 
be observed experimentally through its special spectral 
structure. In fact, the prompt advance in experiments of 
trapping and cooling atoms has provided a platform to 
test our predictions exactly. 

The Hamiltonian for our cold atom system driven by 
two laser beams (see Fig. 1) can be written in the form 
H = P"^ /2M + Hf{r) with the local internal Hamiltonian 



Hf{r)^Y.^,{r)\^) {^\+A\e) (e| 

i=l 

+np (r)|e) {l\+n, (r)|e) {2\ + h.c. 



(1) 



Here, r is the atomic position, A the one-photon de- 
tunning, (r) (Op (r)) the Rabi frequency of the probe 
(control) beam, and Vi (r) the trap potential of the i-th 
inner energy level. We assume that there is no trap for 
the level e, i.e., Ve (r) = 0, and there is the same trap 
for both level 1 and 2, i.e., Vi (r) = V2 (r) = V (r). This 
assumption just ensures the occupance of the dark state. 
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FIG. 1. Three level atoms interacting with two laser beams. 
The probe beam coupling the states |e) and |1) has an or- 
bital angular momentum. The atoms in state |1) and |2) are 
trapped by the potential V, which is a function of the atomic 
position r. 

To use the generalized Born-Oppenheimer approxima- 
tion [8] , we first diagonalize the interaction part Hf (r) of 
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the Hamiltonian and obtain the r-dependent eigenvalues 
of Hf{r): Eo (r) = ^ (r) and 
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E± (r) = -(A ± V'4 + 4 \%\' + K^) + V (r) . (2) 

where A = A (r) = A — (r) is the local one-photon 
detunning. The eigenstate corresponding to Eq (r) is the 
r-dependent dark sate defined as 



|Z?(r)) = l[17p|2)-f7,|l)]. 



(3) 



where VL — ^l{v) = J \VL. 



I rip I . The other two eigen- 



states corresponding to the eigenenergies E± (r) can be 
noted as \B± (r)). Their explicit expressions are not nec- 
essary to the following discussion for the case that the 
atom spatial motion is sufficiently slow so that the inter- 
nal motion is not excited. 

It is well known that the atomic wave function in the 
r-representation can be written as 



(r|vI/)=V.o(r)|i^(r))+ ^ ^u{v)\Bk{v)) 



(4) 



When the energy gaps E± — Eq (r) between the dark 
state \D (r)) and the states \B± (r)) are large enough, 
the Born-Oppenheimer approximation is applicable. Un- 
der this approximation, the atom can be assumed to be 
"kept" in the dark state \D (r)) at every position r and 
the eigen wave function of the total Hamiltonian H can 
be written as (r| = ipo (r) \D (r)) where -00 (r) satisfies 
the eigenequation 



^(P + Aof + y(r) 



V'o (r) = EipQ (r) . 



(5) 



Here, Aq (r) = —i {D (r)| V \D (r)) is the induced gauge 
potential corresponding to the dark state. 

As in Ref. [12], we may assume that the probe beam 
and the control beam have the same frequency and prop- 
agate along the z direction. Then the Rabi frequencies 
rip (r) and flc (r) can be expressed as 



flp — I rip I exp (ik • r + gcj)) 



ric — |ric| exp (ik • r 



(6) 



Here, the real parameters |rip(r)| and |ric(r)| are just 
the slowly varying norms of the fip and Oc and <j) the 
directional angle of the x — y plane. Note that in writing 
down Eq. (6), we have also assumed that the probe beam 
has an orbital angular momentum gcj) with g an integer 
[11]. Then the induced gauge potential Aq of the dark 
state can be written as 



Ao - 



-V0. 



(7) 



It is obvious that the orbital angular momentum num- 
ber g of the probe beam induces an effective magnetic 
monopole of strength g. We thus conclude that the ef- 
fective magnetic charge of this induced monopole can be 
controlled artificially by adjusting the angular momen- 
tum of photons in the probe beam. 

When the norms |Op(r)| and |ric(r)| of the Rabi 
frequencies take the forms |flp|^ = [rid^ = 

^(r — z), the potential Aq has the same form as that 
of the potential created by a monopole and can be ex- 
pressed as 



Ao = V0C 
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cos ( 



2r sin( 



(8) 



Here, r, 9 and (p are just the spherical polar coordinates. 
In this case we have \VLp\ — \Vlc\ — at the origin r = 0. 
This leads to an "accidental degeneracy" E± (0) = Eq (0) 
at the origin where the energy levels cross. Thus the 
Born-Oppenheimer approximation does not work well. 
For this reason, in the following discussion, we only dis- 
cuss the atomic motion around the region far away from 
the origin. 




FIG. 2. The whole space excluding the origin O is divided 
into two regions Ra and /?(,. Ra is the space excluding the 
lower circular cone (the region enveloped by the dashed line) , 
is the space excluding the upper one (the region enveloped 
by the dash dot line). The space outside of the two cones is 
the overlap of Ra and Rh- 

For the original Dirac's monopole, as is well known, 
the induced monopole potential Aq has singularity at 
the string of = 0. But this singularity can be exorcized 
by means of the approach developed by Wu and Yang 
[3,13]. To this end one needs to divide the total real 
space (except the origin) into two regions (see Fig. 2), 
Ra - Q <0 < ^+5 and i?h : f - 5 < 6* < tt, which are 
two overlapping caps and have their well defined local 
coordinates. Then the dark state \D (r)) and the wave 
function %pQ (r) can no longer be defined globally. Instead, 
|D (r)) and ■00 (r) will have different expressions in the 
different regions marked by a and b. For instance the 
dark state can be defined as 

(:,)), = e-^(^=+«^)^(17p|2)-f],|l)) (9) 

in Ra and 
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\D{x)), = er^'^^^{np\2)-n,\i)) (lo) 

in Ri,. The Schroedinger equation (5) can be rewritten 
in different caps as 



^(P + Ao.)Vy(r) 



i^oa (r) = E^boa. (r) (f 1) 



in Ra for a = a,b. Here, (r) and i/'ofc (r) are the 
expressions of ipo (r) in Ra and i?f, respectively and the 
gauge potential Aoa {A.ob) can be expressed as 



Aoa = 



g(-l + cos6') 
2r sin 



g(f + cos6') 
2r sin 9 



60. (12) 



Apparently, Aoa (Aob) is not singular in the region Ra 
(Rb)- In the overlap of Ra and _Rb, we have a connection 
V^ofc (r) = "000 (r) e"*^"^ due to the U{1) — guage. Such 
two local wave functions in two caps with this connection 
in the overlapping region are mathematically called the 
wave section. 

We now consider a simple case that the trap poten- 
tial V (r) is spherically symmetrical. In this case we can 
separate the radical degree of freedom r and the angular 
degrees of freedom 9 and (p in the Schroedinger equa- 
tion (11) by invoking the generalized angular momentum 
operator [13] 



symmetrical, it is convenient to discuss this problem with 
the cylindrical coordinate {p^z,4>). Because of the cylin- 
drical symmetry of V and Aq, the wave function ip^a has 
the factor exp {im(j>) (m — 0, ±1, ±2, ...) which is just the 
generator of the rotation along z axis. Then we have 
ipoa = Tm {p, z) exp {im4>) where Rm (p, z) satisfies the 
radical Schroedinger equation 



1 

2M 



-d''p--dp-d^,+F^{p,z) 
P 



Trr 



VTr^, = ET„ 



Here, the function i^,„ (p, z) is defined as 



Fm (P, z) 



2^p2 ^ ^2p 




(17) 



(18) 



EL 

2r' 



(13) 



It was proved [13] that i/jq (r) = Ri (r) Vg,;,™ can be ex- 
pressed in terms of the monopolc harmonics Y^j^m i^i 4>) 
which is the common " eigensection" of and with 
respect to the eigenvalues I {I + 1) and m. Here, we have 
Z = 1 1 1 , 1 1 1 + 1, ... and m = —I, —1 + 1, I. The radical 
wave function Ri (r) satisfies the equation 



dr{r^dr) , l{l + l)-ig/2f 



2Mr^ 



2Mr^ 



+ V{r)-E 



Ri^O. 



(14) 



In the special case V = 0, when E > 0, the solution of 
Eq. (14) is [14] a Bessel function R = ^= (kr) with 



P = \ll{l + l)-{^ 



i, k = V2ME. (15) 



Next we consider a more interesting case. We assume 
the the trap is a harmonic potential 



V = hlujl {z 
where zq > and p^ — x 



zof + ^Muj'p 



2 „2 



(16) 



;2 +y'^ (see Fig. 3). In this case, 
the atom is confined near a fixed point (0, 0, zq) in the re- 
gion Ra- Therefore, we need only consider the expression 
ipoa of tpo in Ra- Since the trap potential is cylindrically 



FIG. 3. The atoms are trapped in the "discus form" region 
whose center is the point (0, 0,20). In this figure, O is the 
origin of the coordinate. 

If the trap along z axis is strong enough, i.e., lOz is 
large enough, we can make the approximation F^ (p, z) « 
Em (p, zq). It is obvious that F^ (p, zq) can be expanded 
as a Laurent series of p. We also assume that the trap 
localized in x — y plane is also strong enough that we can 
only keep Em {p,zq) up to the term proportional to p^. 
Thus approximately we have 



Tfl 

Em (p, z) « ^ 
P 



16^' 



mg 
2zl' 



(19) 



and we can solve the Eq. (17) to obtain the energy spec- 
trum 



E, 



(2np 



l)w- 



mg 
AMz^ 



(20) 



in terms of the radical quantum number rip and the ver- 
tical one Uz {— 0, 1, 2, ...), where 




9' 



16M2z4 



(21) 



is the modified radical frequency for the two dimensional 
reduced radical oscillator. The corresponding radical fre- 



quency shift uj - 



mg 
' AMz'i 



can be regarded as the first 
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observable effect of the artificial magnetic monopolc. The 
additional term —j^i in the energy spectrum of the 
spatial motion of atom may reflect its effect in realistic 
experiment. The corresponding wave function can also 
be obtained explicitly: 

^m,n„n. ^^^^gi„0^|„|g_M(w+^.)pV2 (22) 

xF {-np, \m\ + l,Mwp2) (^^/Mu^^z^ . 

Here, N^^ = [/MwJ/v^2"-n^!] ^ F is the confluent 
hypergeometric function and iJ„^ the Hermit function. 

Since the above results arc achieved with the gener- 
alized Born-Oppenheimer approximation , we should in- 
vestigate the condition under which this approximation 
is applicable. The adiabatic condition can be obtained 
semiclassically. For simplicity, we only consider the case 
A = 0. In our problem, the suflacient condition of adia- 
batic approximation is 



{D\ V \B±) ■ v| • \E± - Eo\ « 1, 



(23) 



where v is the velocity of the atomic center of mass 
(cm). By straightforward calculation, it can be obtained 
that the upper limit of \{D\ V \B±) ■ v| is ^{\vp\ + \vz\ + 
g\v^\) ^ gv/r. Here, Vp, and are the components 

of the cm velocity in the cylindrical coordinate and v 
the speed rate. On the other hand, the lower limit of 



SnI is 



V 



condition (23) can be rewritten as 



Then the 



l\m«\i^M¥+W+^-E). (24) 

Here we have used the fact that the upper limit of v is 
\J2E/M and the upper limit of V is the atomic energy 
E. 

The adiabatic condition (24) can be satisfied in some 
realistic cases. For instance, we confider the cesium 
atoms trapped around the origin r — 0. We assume 
^ = TT X 10^" Hzm^^/^ and the atomic energy E ^ 10~^^ 
J. Then it is easy to see that when the optical angu- 
lar momentum g ^ 10^, the adiabatic condition is sat- 
isfied when r > 10^^ m. Therefore, if the scale of the 
atomic ensemble is r ^ 10^'^ m, the adiabatic approxi- 
mation holds in almost all the region of the atomic mo- 
tion. Another case discussed above is that the cesium 
atoms are trapped around the point (0, 0, zq) by the har- 
monic potentials. We assume zq ~ 10~^ m, 5 ~ 10'' and 
the frequencies of the trap potentials are Uz ^ 10^ Hz 
and w ~ 10^ Hz. Then the adiabatic condition can be 
met again if ^ and E have the same values as mentioned 
above. Therefore, the induced change of zero point en- 
ergy —mg/ (4Mzq) is about 10^ Hz. The frequency shift 
ui — IV caused by the monopole potential will have the 
same order with co. 



Before concluding this paper we make a remark on the 
results that we have obtained. As we have shown, when 
the certain conditions = ^(r + z), \0,cf = — z) 
is satisfied, the singularity of the monopolc potential can 
be eliminated by dividing the space into two caps. To 
see the universality in choosing the profile of light beam 
we make the following change: jfJcl^ = C [(2?? — 1) r — z] 
where 77 > 1. In this case, the induced gauge potential 
will be Ao = g {1 + cos6) e^/{2r]rsm9). Thus the sin- 
gularity of Aq can NOT be eliminated with the above 
method when g/rj is not an integer. This is because in 
this case we will have V06 (r) = V'Oa (r) e"'^'^/'' and if 
g/rj is not an integer, e~^^'^/^ is not a single valued func- 
tion. We notice that a similar situation appears when 
the Dirac's quantization condition is not satisfied [13]. 
This concluding remark shows the possibility to test the 
subtle effect of magnetic monopole with manipulation of 
ultra cold atoms. 
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